(Preprint) A AS 1 1 -XXX 


MULTI-MANEUVER CLOHESSY-WILTSHIRE TARGETING 

David P. Dannemiller 


Orbital rendezvous involves execution of a sequence of maneuvers by a chaser 
vehicle to bring the chaser to a desired state relative to a target vehicle while 
meeting intermediate and final relative constraints. Intermediate and final relative 
constraints are necessary to meet a multitude of requirements such as to control 
approach direction, ensure relative position is adequate for operation of space-to- 
space communication systems and relative sensors, provide fail-safe trajectory fea- 
tures, and provide contingency hold points. The effect of maneuvers on constraints 
is often coupled, so the maneuvers must be solved for as a set. For example, ma- 
neuvers that affect orbital energy change both the chaser’s height and downrange 
position relative to the target vehicle. Rendezvous designers use experience and 
rules-of-thumb to design a sequence of maneuvers and constraints. 

A non-iterative method is presented for targeting a rendezvous scenario that in- 
cludes a sequence of maneuvers and relative constraints. This method is referred to 
as Multi-Maneuver Clohessy- Wiltshire Targeting (MM_CW_TGT). When a single 
maneuver is targeted to a single relative position, the classic CW targeting solution 
is obtained. 

The MM_CW_TGT method involves manipulation of the CW state transition 
matrix to form a linear system. As a starting point for forming the algorithm, 
the effects of a series of impulsive maneuvers on the state are derived. Simple 
and moderately complex examples are used to demonstrate the pattern of the re- 
sulting linear system. The general form of the pattern results in an algorithm for 
formation of the linear system. The resulting linear system relates the effect of ma- 
neuver components and initial conditions on relative constraints specified by the 
rendezvous designer. Solution of the linear system includes the straight-forward 
inverse of a square matrix. Inversion of the square matrix is assured if the de- 
signer poses a controllable scenario - a scenario where the the constraints can be 
met by the sequence of maneuvers. Matrices in the linear system are dependent 
on selection of maneuvers and constraints by the designer, but the matrices are 
independent of the chaser’s initial conditions. For scenarios where the sequence 
of maneuvers and constraints are fixed, the linear system can be formed and the 
square matrix inverted prior to real-time operations. 

Example solutions are presented for several rendezvous scenarios to illustrate 
the utility of the method. The MM_CW_TGT method has been used during the 
preliminary design of rendezvous scenarios and is expected to be useful for itera- 
tive methods in the generation of an initial guess and corrections. 


’Aerospace Technologist. Flight Dynamics Division, DM34, NASA-Johnson Space Center, 2101 NASA Parkway, Hous- 
ton, Texas, 77058, David.P.Dannemiller@nasa.gov. 
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A non-iterative method is presented for targeting a rendezvous scenario that in- 
cludes a sequence of maneuvers and relative constraints. This method is referred to 
as Multi-Maneuver Clohessy- Wiltshire Targeting (MM_CW_TGT). When a single 
maneuver is targeted to a single relative position, the classic CW targeting solution 
is obtained. The MM_CW_TGT method involves manipulation of the CW state 
transition matrix to form a linear system. Solution of the linear system includes 
the straight-forward inverse of a square matrix. Example solutions are presented 
for several rendezvous scenarios to illustrate the utility of the method. 

INTRODUCTION 

Orbital rendezvous involves execution of a sequence of maneuvers by a chaser vehicle to bring 
the chaser to a desired state relative to a target vehicle while meeting intermediate and final relative 
constraints. Intermediate and final relative constraints are necessary to meet a multitude of require- 
ments such as to control approach direction, ensure relative position is adequate for operation of 
space-to-space communication systems and relative sensors, provide fail-safe trajectory features, 
and provide contingency hold points. The effect of maneuvers on constraints is often coupled, so 
the maneuvers must be solved for as a set. For example, maneuvers that affect orbital energy change 
both the chaser’s height and downrange position relative to the target vehicle. Rendezvous designers 
use experience and rules-of-thumb to design a sequence of maneuvers and constraints. 

A non-iterative method is presented for targeting a rendezvous scenario that includes a sequence 
of maneuvers and relative constraints. This method is referred to as Multi-Maneuver Clohessy- 
Wiltshire Targeting (MM_CW_TGT). When a single maneuver is targeted to a single relative posi- 
tion, the classic CW targeting solution is obtained. 1 

As a starting point for forming the MM_CW_TGT algorithm, the effects of a series of impulsive 
maneuvers on the state arc derived. The MM_CW_TGT method involves manipulation of the CW 
state transition matrix to form a linear system. Simple and moderately complex examples are used 
to demonstrate the pattern of the resulting linear system. The general form of the pattern leads 
to an algorithm for formation of the linear system. The resulting linear system relates the effect 
of maneuver components and initial conditions on relative constraints specified by the rendezvous 
designer. Solution of the linear system includes the straight-forward inverse of a square matrix. 
Inversion of the square matrix is assured if the designer poses a controllable scenario - a scenario 
where the constraints can be met by the sequence of maneuvers. 
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ton, Texas, 77058, David.P.Dannemiller@nasa.gov. 


1 



A numerical example is presented to illustrate the utility of the method and to provide comparison 
data. 

NOTATION 

Vectors arc denoted by an over arrow (e.g. r). Subscripts denote time or event information. For 
example, ro is the value of r at time to. 




Figure 1. Local Vertical Curvilinear Coordinate System 

State vectors are in Local Vertical Curvilinear (LVC) coordinates (axes orientation is shown in 
Figure 1) and arc represented as: 


s = 


r 

v 


where r is the position vector and v is the velocity vector. 



X 

,v = 

X 

r = 

y 

y 


z 


z 


The 6x6 Clohessy-Wiltshire (CW) transition matrix is denoted by 6'” where the subscript denotes 
the initial time and superscript denotes the final time, t rn and t n respectively in this case. 

Impulsive maneuvers arc instantaneous changes in velocity at a specified time. To remove ambi- 
guity, it is necessary to specify pre or post-maneuver when referring to velocity and state vectors. 
Subscripts - and + arc used to refer to pre and post-maneuver quantities respectively. 


AF n = v n+ - v n - 


A y n 

A z n 


While not strictly required, subscripts - and + arc also applied to times to differentiate an interest 
in pre or post-maneuver quantities (e.g. i n _, t n +). 

EFFECT OF A SERIES OF IMPULSIVE MANEUVERS 

As a starting point for forming the MM_CW_TGT algorithm, the effects of a series of impulsive 
maneuvers on the state arc derived. 
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Given an initial state at time to- and impulsive maneuvers AFo, Av\, and AF 2 at times to, 
1 1 , and t 2 , determine the state at time t 2+ . 


s*o+ 

si- 

S1+ 

S2+ 


, 0 

s °- + . ^ 
[ZX^oJ 

Cl • so+ 

Cq ■ so- + Cl ■ 


0 

Ahj 


0 

AF 0 


0 

An 0 


si- + 

Cl ■ so_ + Cl ■ 

Cl • si+ + 

Cl ■ (cl- so- + Cl 


0 

Av 2 


Cq ■ So- + Cn 


0 

_Av 0 _ 


+ 


0 

Avi 


0 

Av 0 

+ cl 


+ 


0 

Av\ 

0 ' 
Av\ 


+ 


0 

Av 2 


+ 


0 

Av 2 


In general there will be impulsive maneuvers at times to, t\, t 2 , . . . t n . The general form for the 
effect of this series of n impulsive maneuvers is obtained by extrapolating the pattern above. 


S n + = C% ■ So- + 


1 oi 
1 1 

+ c? • 

0 ' 
Av\ 

+ C%- 

0 ' 
_Av 2 _ 


.. + 

't> 

^ 01 
L 1 


(1) 


The MM_CW_TGT method involves manipulation of the CW state transition matrix to form a 
lineal - system. Two examples are used to demonstrate the pattern of the resulting linear system, a 
simple targeting example and a more complex targeting example. 


SIMPLE TARGETING EXAMPLE 

Assume a pair of horizontal maneuvers transfer a chaser vehicle from a stable location on the x 
axis at s*o_ to another stable location Sf + in an integral number of orbits. Figure 2 shows the desired 
trajectory. 

For this example, Eq. (1) takes the following form. 


»/+ 


C f o • so- + Cl 


0 ' 
Avo 


+ 


0 

A ^>. 


(2) 


Restructure this equation in 2 steps. The first step takes advantage of knowledge of the selected 
maneuver components. Since Avq and A vj are specified as horizontal, they take the form 


Av 0 = 

1 

O 

• B o 

<1 

1 

ii 

< 

1 


L o J 


L o J 
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Figure 2. Simple Targeting Example 


Substitute Eq. (3) into Eq. (2). 
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* 
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* 
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0 
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0 
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0 
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* 
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One of the CW transition matricies is expanded. Most of the terms are denoting that their 
specific values do not matter. Because of the form of Avq, only one column is relevant. The identity 


matrix is added in front of the Av f term for symmetry. Simplifying, 


«/+ 


«/+ 


Cq ■ -so- + 


C'o ‘ s o- + 


^(1,4) 
CU 2,4) 
Co (3,4) 
^(4,4) 
(5)4) 
^(6,4) 



O' 


0 

• Ax 0 + 

0 

1 


0 


_ 0 _ 


Co (1,4) 0 
CU 2,4) 0 
(3)4) 0 
^(4,4) 1 
C'o (5)4) 0 
.C'o (6,4) 0. 


Ax 0 

A®/. 


• Axy 


(4) 


The second restructuring step takes advantage of the selected constraints. For this example the 
final state is specified as a stable location on the x axis (i.e. if = 0). This implies that Sf + takes 
the form 


x f 

* 



* 

* 


(5) 


Note that yj, Zf, iff, and zf are not explicitly constrained. Given the example setup (i.e. initial 
state is stable on the x axis, and transfer time is an integral number of orbits), by definition zj and 
Zf are both zero. Since no maneuver components are specified in the y direction, neither yf nor ij f 
can be used as constraints. 

Substitute Eq. (5) into Eq. (4), and then remove all but the Xf and if rows. 


V 


c 0 f ( 1,1) 

cl( 1 , 2 ) 

Co (1, 3) 

C 0 / (l,4) 

Co (1, 5) 

Co (1,6)" 

0 


pi (4, 1) 

C$( 4,2) 

Cq (4, 3) 

Cq (4, 4) 

Cq (4, 5) 

C 0 / (4,6)_ 


"Cq (1,4) o' 


Ax 0 

_Co (4,4) 1. 


A*/. 


Eq. (6) is a system of linear equations with 2 equations and 2 unknowns - the desired maneuver 
components Axo and A if. Use shorthand notation to compact Eq. (6). 

S constraint — ^4 • S()— T B • Av (7) 

The n-vector s cons t r aint contains the selected constraints, where n is the number of constraints. 
The nx6 matrix A is the effect of sq- on the constraints. The n-vector Av contains the selected 
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maneuver components. The nxn matrix B is the effect of Av on the constraints. Rearrange Eq. (6) 
to solve for Av. 


Av — B • ( Sconstraint A • SO— ) (8) 

Applying the maneuver components Axo and A x/, at the specified times, will satisfy the con- 
straints. 

To summarize, a sequence of maneuvers (times and components) and constraints (times and 
components) were designed. A system of linear equations was structured in 2 steps. The first 
step eliminated terms associated with zero maneuver components. The second step formed the 
Sconstraint. vector, and the A and B matrices by keeping only the terms associated with constraints. 
Since the number of maneuver components is equal to the number of constraints, the resulting B 
matrix is square. Inversion of the B matrix is assured if the designer poses a controllable scenario - 
a scenario where the the constraints can be met by the sequence of maneuvers. 

MODERATELY COMPLEX TARGETING EXAMPLE 

Using the multi-maneuver CW technique to solve the simple example in the last section is 
overkill, but provides initial insight into the pattern of the resulting linear - system. More insight 
is provided in this section by a moderately complex example. 



Figure 3. Moderately Complex Targeting Example 

For this scenario, a sequence of maneuvers NC, NH, TI (Transition Initiate), TF (Transition Final) 
transfer a chaser vehicle from an initial state on the x axis at sq-, to a desired TI position at sti+, 
to the final state at stf+ ■ Figure 3 shows the desired trajectory. 

The burn names used here have Gemini, Apollo, and Shuttle heritage. “N” in the burn names 
was originally (circa 1964) a Mission Control Docking Initiation (DKI) rendezvous burn targeting 
program counter variable. “N” stood for the N th crossing of the chaser line of apsides where the burn 
was performed (such as 1 for first apogee, 1.5 for first perigee, 5 for fifth apogee, etc.). The burn 
type was added using subscripts as lc (Catch-up or phasing), 2.5h (Height), 3pc (Plane Change), 
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4sr (Slow Rate or coelliptic), and 4.5cc (Corrective Combination). Eventually the apse count and 
subscripts notation were dropped. The resulting burn names became NC, NH, NPC, NSR, and NCC. 

We specify horizontal maneuvers at NC and NH, and 2-axis maneuvers at TI and TF. Typically, 
TI and TF include maneuver components in the y direction. They arc not included in this example 
in order to keep the problem tractable. The maneuvers take the following form. 


Avnc = 

Ax N c 

0 

, Av NH = 

Ax nh 

0 

II 

<1 

A XTI 
0 

, Av tf — 

Axtf 

0 


0 


0 


.A ZTI_ 


Aztf 


We have constraints at TI and TF. 


~X T { 


xtf 

* 


* 

zti 

* 

, STF+ = 

ztf 
xtf- i- 

* 


* 

* 


_ZTF+_ 


Notice that the number of maneuver components and the number of constraints arc both 6, which 
will result in a square B matrix. Also notice that the scenario is controllable; NC/NH arc capable 
of maneuvering the chaser to the specified TI position xti and zti '■ TI/TF are capable of achieving 
the desired final conditions. 

Since constraints arc specified at 2 times, the effects of maneuvers at both TI and TF must be 
formulated. First TI. 


STI+ — Cq 1 • So- + Cj/tQ 


0 

A v NC 


, r TI 

+ L 'NH 


0 

A v N h 


+ 


0 

Avj ’i 


Next TF. 


? — n TF 

STF+ — 


so - + C t n F c 


0 

_Av nc _ 

i r TF ■ 
T- f 'NH 

0 

_Av nh _ 

i r TF ■ 

T- O'T/ 

> 

1 1 

+ 

0 

A vtf 


The equations are restructured in 3 steps. The first step takes advantage of knowledge of the 
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selected maneuver components. 


-* s~iT T —> 

STI+ = C 0 ■ Sq- + 
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0 
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Cyf (5, 6) 
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'0 

O' 



0 

0 


Ax T / 

+ 

0 

0 
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1 

0 



0 
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_0 

1 


A xtf 
A.Ztf 


A XT/ 
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The second step in restructuring takes advantage of the selected constraints. Eliminate all rows 
that arc not constraints. 


X T / 


_2T/ 



<#'(M) 

.<#'( 3,1) 


<#'( 1 , 2 ) 

C 0 ^(3,2) 


<#'(1,3) C#'(l,4) 
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<#'( 1,6/ 
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For the final step in restructuring, stack the 2 equations into a single linear system. 
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CVh( 4,4) 

C£f( 4,4) 

C£f( 4,6) 

1 

0 


A xtf 

4) 

Cnh(6, 4) 
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• so- 


(9) 


Notice that Eq. (9) has the same general form as Eq. (7). We solve for Av in the same manner as 
Eq. (8). 

GENERAL PATTERN 

We use Eq. (9) to find the general pattern of the Multi-Maneuver CW Targeting algorithm. 


1 . The selected maneuver components and constraints are used to form a system of linear equa- 
tions of the form Sconstraint = A ■ s 0 - + B ■ Av. 

2. The number of maneuver components must be equal to the number of constraints, which is 
denoted as n. 

3- s cons traint is an n-vector that contains the selected constraints. 

4. Av is an n-vector that contains the selected maneuver components. 

5. A is an nx6 matrix that contains the effects of so- on die constraints. Each row of A contains 
a row from the CW transition matrix, the specific row being that associated with the selected 
constraint, and the transfer time equal to t cons t ra int ~ to- 

6. B is an nxn matrix that contains the effects of the maneuver components on the constraints. 
Each row of B contains terms from the CW transition matrix, the specific row being that 
associated with the constraint, the specific column being that associated with the maneuver 
component, and the transfer time equal to t cons t ra int — tmaneuver ■ If the constraint time is 
after the maneuver time, the term is 0 (later maneuvers cannot effect earlier constraints). 

7. B must be invertible. This is assured if the designer poses a controllable scenario - a scenario 
where the the constraints can be met by the sequence of maneuvers. 

8. The linear system is solved for the unknown maneuver components. 

Av — B ( Sconstraint A • SQ— ) 
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ACCELERATION CONSTRAINT 


Many trajectory scenarios can be solved with the position and velocity constraints discussed 
above. Two other useful constraint types arc derived in this section and next. 

Co-elliptic segments are useful in many rendezvous scenarios. Targeting co-elliptic conditions 
with only position and velocity constraints is possible, but awkward. When co-elliptic, z = const, 
z = 0, and z = 0, so the addition of an acceleration constraint makes co-elliptic targeting straight 
forward. 

The CW differential equations for coasting flight arc: 


x = 2 • uj ■ z 
V = -u 1 ■ y 
z = 2> ■ uj 2 ■ z — 2 • ui ■ x 


Put these equations in matrix form. 


X 


'0 

0 

0 

0 

0 

2 • uj 


r 

y 

— 

0 

-cu 2 

0 

0 

0 

0 


z 


0 

0 

3 • tc 2 

—2 • uj 

0 

0 


V 


'a = PVA ■ 


r 

v 


( 10 ) 


We refer to the 3x6 matrix in Eq. (10) as PVA since it converts a position/velocity vector to an 
acceleration vector. Notice that PVA is a function only of uj; it is time independent. Pre-multiply 
Eq. (1) by PVA. 


fln+l — PVA • S n -|- 


PVA ■ C% • so- + PVA ■ C'S • 


0 

Av 0 


+ . . . + PVA ■ 


0 

Av n 


( 11 ) 


For scenarios where an acceleration constraint is desired, form the constraint using Eq. (11) and 
adjoin it to the linear system. 

ELEVATION ANGLE CONSTRAINT 

Another useful constraint is elevation angle. Figure 4 shows the convention for the in-plane 
elevation angle. 

Specifying an elevation angle constraint imposes a constraint on the relationship between x and 
If 6 is the desired elevation angle, then 


Form the following. 


, n ze 

tan 6 = 

x e 


xg = f x ■ so + 9x ■ Av 
zg = fz ■ s 0 + g z ■ Av 


( 12 ) 


(13) 

(14) 
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(-)-• 


N 0 - 


(+)- 


(+) 0 (-) 
x 


Figure 4. Elevation Angle Convention 



Here f x and g x arc row vectors that represent the effect of sq and An on x$. Likewise for f z and 
g z on zq. Substitute Eqs. (13) and (14) into Eq. (12) and restructure. 


0 = (j z + tan# • 


■ so + (g z + tan 0 • g x ) ■ Av 


(15) 


Note that Eq. (15) has a form similar to previous constraints. For scenarios where an elevation 
angle constraint is desired, form the constraint using Eq. (15) and adjoin it to the linear system. 


NC, NH, NSR, TPI, TPF NUMERICAL EXAMPLE 

This section provides a numerical example to illustrate the utility of the method and to provide 
comparison data. Figures 5 and 6 illustrate the desired relative motion. 

Initial conditions for the scenario are near the x axis with some out-of-plane position. 


to = 1000 sec, 


So = 


- 1 , 000 , 000 ' 

2000 

0 

0 

0 

0 


ft 


The target orbit period is 90 minutes. 


2 • 7 r rad 
5400 sec 
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Figure 5. Numerical Example 


Event times are as follows. 


to = 1000 sec 
tNC = to 


IT 


tNH = tNC + 2.5 • — 

W 

tNSR = t-NH + .4 • — * 

W 
7 r 

trpi = tNSR + -5 • — 
w 

130 vr 

t T PF = tTPI + NFF ' ~ 

360 W 


( i.e . 2.5 orbits from NC ) 


(i.e. 130 degree transfer from TP I) 


Maneuvers NC and NH are horizontal. NSR has x and z maneuver components. TPI and TPF 
are 3-axis - x, y, and 2 maneuver components. 


Av nc = 


A x N c 


Ax nh 


Axnsr 

0 

, Av N h = 

0 

, A vnsr = 

0 

0 


0 


_Aznsr_ 



A XTPI 


A xtpf 

Avtpi = 

AijTPi 

, Avtpf = 

AijTPF 


_Aztpi_ 


Aztpf 


Constraints are imposed at 3 event times - post-NSR co-elliptic, TPI elevation angle, and post- 
TPF state. 


ZNSR = 10,000 ft, ZNSR+ = 0, ZNSR+ = 0 
Otpi = -5 rad 
STPF+ = 0 

‘Using a transfer of .5 orbits instead of .4 would result in lower Av. This example uses .4 to better illustrate the NSR 
maneuver. 
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Figure 6. Numerical Example (cont.) 


We have everything required to form the linear system, Eq. (7). Doing so results in the following 
quantities. 


S constraint — 


ZNSR 

ZNSR+ 

ZNSR+ 

row associated with dp pi 
XTPF+ 

VTPF+ 

ZTPF+ 

XTPF+ 

VTPF+ 

ZTPF+ 


Av = 


' Axatc ' 

Ax nh 

A x NS r 
Aznsr 
A xppi 

A i)TPI 

Aztpi 

Axt p p 
A ijTPF 

.A ZppF . 
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Solve this using Eq. (8). 


A VNC 


A V NH 


Avnsr 


A vt pi 


L \vtPF 


-18.96 

0 

0 


- 1.21 

0 

0 

14.36' 

0 

20 . 86 _ 

' 3.05 ‘ 
- 0.21 
-2.09 


sec 

ft_ 

sec 

A 

sec 

A 

sec 


2.76 

-2.46 

3.33 


A 

sec 


Applying these A vs at the appropriate times results in the relative motion trajectory shown in 
Figures 5 and 6. 
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CONCLUSION 

When the accuracy associated with CW equations is adequate, the MM_CW_TGT method is 
useful as is. The non-iterative nature of the method results in fast and robust results. The author has 
used MM_CW_TGT for preliminary design of rendezvous scenarios. 

Although not yet implemented, the author expects MM_CW_TGT to be useful as a paid of iterative 
methods. MM_CW_TGT can generate initial guesses for maneuvers, and can calculate corrections 
between iterations. 

ACKNOWLEDGMENT 

The author thanks Aaron Brown for suggestions that improved the readability of the paper. 

REFERENCES 

[1] W. H. Clohessy and R. S. Wiltshire, “Terminal Guidance System for Satellite Rendezvous,” Journal of 
the Aerospace Sciences , Vol. 27, No. 9, 1960, pp. 653-659. 


15 



Multi-Maneuver 

Clohessy-Wiltshire 

Targeting 


David P. Dannemiller 
NASA - JSC 
8 / 4/2011 

(Astrodynamics Specialist Conference) 


Introduction 


• Orbital rendezvous 

- Sequence of maneuvers 

- Bring the chaser to a desired state relative to a target 

- Meet intermediate and final relative constraints 

• Effect of maneuvers on constraints is often coupled, so the 
maneuvers must be solved for as a set 

• Multi-Maneuver Clohessy-Wiltshire Targeting 
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Notation 


• Local Vertical Curvilinear (LVC) 






- 


X 


X 

• State vector 

S = 

r 

V 

f = 

y 

z 

,V = 

y 

z 


• Clohessy-Wiltshire (CW) transition matrix c ™ 

- Subscript - initial time 

- Superscript - final time 


• Impulsive maneuvers 


- (-) pre-maneuver 

- (+) post-maneuver 







^■X n 

^■ijn 

^■Z n 


Flight 

Dynamics 
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Effect of a Series of 
Impulsive Maneuvers 


Two maneuvers 


So- + 

ll 


S 0 + — 

Si- = Cq • s 0 + 
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Cl ■ s 0 - + Cl 


si+ — si_ + 
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A^o 


Cl ■ sq - + Cq 
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A^o 
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A^i 


General 
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- CZ • so- + c. 
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0 


o 

+ CT • 

o n 

+ C? • 

o 

+ . . . + 

o 

AiTo 

i ^ 1 

AiTi 

1 z 

A{T 2 


Av n 


-Li 
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Simple 

Targeting Example 


• Maneuvers 


<1 

1 

o 

o 

<1 

i 

,Av f = 

'A xf 
0 


o 


O 


• Constraints 

~ x f 



Dynamics 
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Simple 

Targeting Example (cont) 
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Simple 

Targeting Example (cont) 


S f+ = 


Cq ■ SQ— + 



Ax 0 

Ax f _ 


cl( 1 , 1 ) c£( 1 , 2 ) Cq ( 1 , 3 ) Cq ( 1 , 4 ) Cq ( 1 , 5 ) Cq ( 1 , 6 )" 

^o(4>3) ^o(4>4) ^o(4>5) ^o(4>6) 
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• So- 
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1 o 
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•H 
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o° 
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•H 
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^ constraint — ^ ' ^ 0 — ~t~ B • Av 
Av = B • (s constraint A ’ SQ — ) 
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Moderately Complex 
Targeting Example 


• Maneuvers 


Avnc = 

Ax nc 

0 

, = 

Ax nh 

0 

, = 

i 

> 

o 

1 

, /S.vj'p = 

Axtf 

0 


0 


0 


A zti 


Aztf 


• Constraints 


STI+ = 


XT I 


XTF 

* 


* 

Z TI 

* 

> s T f+ — 

ZTF 

xtf+ 

* 


* 

* 


_ztf+_ 
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Moderately Complex 
Targeting Example (cont) 


sti+ = Cq 1 ■ sq_ + C 
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Moderately Complex 
Targeting Example (cont) 
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General Pattern 


1 . The selected maneuver components and constraints are used to form a system of linear equa- 
tions of the form s constra int = A ■ s 0 _ + B ■ Av. 

2. The number of maneuver components must be equal to the number of constraints, which is 
denoted as n. 

3- s cons traint is an n-vector that contains the selected constraints. 

4. Av is an n-vector that contains the selected maneuver components. 

5. A is an nx6 matrix that contains the effects of so- on the constraints. Each row of A contains 
a row from the CW transition matrix, the specific row being that associated with the selected 
constraint, and the transfer time equal to t cons t ra int — to- 

6. B is an nxn matrix that contains the effects of the maneuver components on the constraints. 
Each row of B contains terms from the CW transition matrix, the specific row being that 
associated with the constraint, the specific column being that associated with the maneuver 
component, and the transfer time equal to t cons t ra int ~ tmaneuver • If the constraint time is 
after the maneuver time, the term is 0 (later maneuvers cannot effect earlier constraints). 

7. B must be invertible. This is assured if the designer poses a controllable scenario - a scenario 
where the the constraints can be met by the sequence of maneuvers. 

8. The linear system is solved for the unknown maneuver components. 

Av — B • (s constraint A ■ Sq— ) 




Other Constraints 


• Acceleration Constraint 

- Useful to specify co-elliptic constraint 

• Elevation Angle Constraint 

- Often used in terminal phase 


Dynamics 


-Li 
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NC, NH, NSR, TPI, TPF 
Numerical Example 


Maneuvers 


Avnc = 


Ax N c 


Axnh 


0 

, AtTtv// = 

0 

> Av N sr = 

0 


0 



Axnsr 

0 

Aznsr 



A xtpi 


Axrpf 

A vjppj = 

A yrpi 

, Avtpf = 

A vtpf 


Aztpi 


A ZfpPF 


0 -1 -2 -3 


-5 -6 -7 -8 -9 


-10 
x 10 5 


Constraints 

ZNSR = 10,000 ft, ZNSR+ = 0, ZNSR+ 
Otpi = -5 rad 
stpf+ = 0 

Numerical solution in paper 


= 0 
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Conclusion 


• Non-iterative nature of the method results in fast and robust results 

• Method has been used for preliminary design of rendezvous 
scenarios 

• Method expected to be useful as a part of iterative methods 

- Generate initial guesses for maneuvers 

- Calculate corrections between iterations 
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